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The nonequilibrium phase transition of an open Takayama-Lin Liu-Maki chain coupled 
i— i . with two reservoirs is investigated by combining a mean-field approximation and a formula 

' characterizing nonequilibrium steady states, which is obtained from the algebraic field theo- 

O ' retical approach to nonequilibrium statistical mechanics. When the bias voltage is chosen as 

a control parameter, the phase transition between ordered and normal phases is found to be 
of first or second order. The current- voltage characteristics are S-shaped in some parameter 
I ' region. In contrast, when the current is chosen as a control parameter, all the nontrivial solu- 

\ tions of the self-consistent equation are found to be stable. In this case, the phase transition 

. between the ordered and normal phases is always of second order, and negative differential 

^ ' conductivity appears at low temperature. 

£3 . §1. Introduction 

■ 

One of the main interests in mesoscopic physics is a full understanding of the 
nonequilibrium properties of small quantum systems coupled with large reservoirs. 
' When the reservoirs have different temperatures and/or chemical potentials, the 

whole system can be in a steady state with constant particle and energy flows, i.e. 
a nonequilibrium steady state (NESS). One of promising approaches in dealing with 
f*"*- I such systems is based on algebraic quantum field theory.®'®'® For example, with 

the aid of the algebraic approach, Pusz and Woronowicz rigorously derived Carnot's 
formula^® and Oj ima et al. proved the positivity of the (relative) entropy produc- 
tion rate,® both for systems of infinite degrees of freedom. In addition, for the XY 
• model, a NESS was rigorously constructed.® Recently, starting from Ruelle's work 

<3\ . on scattering-theoretical characterizations of NESS® and Jaksic-Pillet's investiga- 

tion into entropy production,® the algebraic approach to NESS has been extensively 
developed [see Refs 15]) .[TO ]) .[TT ]) . [T3j) and references therein]. Currently, linear re- 
sponse theories^® thermodynamic properties^® '^® and the Landauer-Biittiker for- 
' mulcP®'^® are investigated, in addition to various other aspects P® 

Hereafter, we consider a quantum system coupled with two free fermionic reser- 
voirs L and R described by annihilation operators ak ff and &k<r, respectively, where 
k refers to the wave number and a to the spin. In this case, if the reservoirs are 
initially set to be in different equilibria, the whole system is shown to approach a 
NESS in the long time limit provided that the incoming fields of and (3^ 
of frko- are compieteP®'^® The NESS so obtained can be characterized as a state 
satisfying Wick's theorem with respect to a^ a and (5^ a and having the two-point 
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functions: 

( a L a kV')oo = fli^kL^aa' £(k-k') , (/^/W}^ = fR{tw kR )5 a(T i5{\s.-\s!) , (1-1) 

where (• • • )oo stands for the average with respect to the NESS, hujku is the single- 
particle energy of wave number k, f u (x) = l/^ x ~^^ T,J + 1) is the Fermi distri- 
bution function, T v is the initial temperature and ji v is the initial chemical poten- 
tial of the reservoir v = L,R. Formally, this can be understood as follows:*) Let 
po be an initial density matrix, where the two reservoirs are in distinct equilib- 
ria, and poo be that of the NESS, then, \\m t ^+oo e~ tHt f h poe lHt / h = and, e.g., 
r ^{ a k t 7°k'<T'Po} = /L(^/cL)<5o-o-'<5(k — k'). As the incoming field a\^ a is given by 
lim^ +00 e^khi-t) e iH{-t)/h a ^ e -tH(-t)/h _ a ^ one OD tains the desired relation: 

/&(fookL)<W*(k - k') = Tr^ava'PoW^-"*'^ 

= Tr{{e i f(- t )a kCT e- l f(~ t ) e ^(- i )}te i f(-')a kV e- i f(- t )e^'^- i )e- i fVoe iff } 

Tr {«L a k' CT '/Ooo} = (a kCT a k / CT /)oo (as t -> +00) . (1-2) 

As noted by Blanter and Buttiker^S 1 [cf. Eqs.(29) and (36) in their paper], the NESS 
characterization can be a starting point for the Landauer-Biittiker approach to 
transport properties of mesoscopic circuits. Note also that can be applied even 
to systems with interacting fermions if the incoming fields are complete.^ 1 Indeed, 
Katsura successfully derived a NESS from for a solvable model of the Kondo 
effectP 

Since eq. (|l-ip fully characterizes NESSes of noninteracting fermions, it is natu- 
ral to consider a mean-field approximation based on (ll-ip for a NESS of interacting 
fermions. Based on this view, we investigated a NESS of an Aharonov-Bohm ring 
with a quantum dot within a mean-field approximation and obtained a differential 
conductivity consistent with numerical renormalization group analysis and experi- 
mentsPD Here, with the aid of a similar nonequilibrium mean-field approximation, 
we study the nonequilibrium phase transition in the Takayama-Lin Liu-Maki chairP^J 
(TLM chain) embedded between two infinitely extended reservoirs. The TLM chain 
is a continuum limit of a lattice model (the SSH lattice) for polyacetylene proposed 
by Su, Schrieffer and Heeger^J and describes the charge density wave commensurate 
with the lattice. 

We employ the TLM chain as a representative example of systems with phase 
transitions. However, since the mean field approximations of the TLM chain, a su- 
perconductors, a ID extended Hubbard lattice of spinless fermions and the Jordan- 
Wigner-transformed XXZ model are equivalent, the present analysis would provide 
some insight into nonequilibrium properties of various interacting systems, such as 
the current-induced suppression of the charge order observed in some #-type BEDT- 
TTF organic conductors^ 1 J25},[26},[271i an( j negative differential conductivity re- 
cently reported for the XXZ modeP^ and some strongly correlated systems 

The rest of this paper is arranged as follows. In Sec. 2, we introduce a finite 
TLM chain coupled with two infinitely extended reservoirs. In Sec. 3, a mean-field 

The very proof of the existence of the limits requires rigorous and careful arguments.^ 
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approximation based on (|lTj) is discussed. The averaged lattice distortion serves 
as an order parameter, and its self-consistent equation is obtained by averaging the 
equation of motion of the lattice distortion with respect to a nonequilibrium steady 
state. In Sec. 4, the self-consistent equation, current and stability conditions are 
explicitly derived in case where the TLM chain is long enough and the order pa- 
rameter is spatially uniform. In Sec. 5, possible phases are discussed in detail when 
the chain-reservoir couplings are symmetric. When the bias voltage is chosen as one 
of the control parameters, the phase transition between ordered and normal phases 
could be of first or second order depending on the bias voltage and temperature. At 
low temperature, the current-voltage characteristics are S-shaped. For some bias- 
voltages, the temperature dependence of the order parameter is found to be similar 
to that for the nonequilibrium superconducting phase induced by excess quasiparti- 
clesP^'''® In contrast, when the current is chosen as one of the control parameters, 
the self-consistent equation has a unique stable solution and the phase transition be- 
tween the ordered and normal phases is always of second order. Negative differential 
conductivity appears when the temperature is lower than a certain threshold value. 
In Sec. 6, after a summary of the paper is given, the self-consistent equation for 
the open TLM chain is compared with that for the nonequilibrium superconductor 
obtained by Scalapino et alW^ Then, on the basis of the similarity of mean-field 
approximations of the open TLM chain and an open ID extended Hubbard lattice, 
the experimental results for some #-type BEDT-TTF organic conductors^ 
as well as the negative differential conductivity found in an open XXZ modePSl 1 are 
qualitatively discussed within the scope of the present analysis. In Appendix A, 
an open TLM chain is derived from an open SSH lattice. In Appendix B, normal- 
mode operators are explicitly given. In Appendix C, we discuss the relationship 
between the average chemical potential and the Coulomb energy. In Appendix D, 
Green functions necessary for deriving the self-consistent equation are provided. In 
Appendix E, the stability of the nontrivial phases at zero temperature is discussed. 
In Appendix F, the Ginzburg-Landau expansion for the self-consistent equation is 
given. 



§2. Open TLM Model 



The system in question consists of a finite TLM chain and two free electron 
reservoirs. In terms of the quantized local lattice distortion A(x) and the two- 
component electron field <2v(x) 



d a (x) 
e a (x) 



the Hamiltonian of the TLM chain is given b 



„ J o 







-ihva v — — h A(x)a x 
y dx 



WJx) 



(2-1) 



1 

2irhvX 



dx 



A{x) 2 + -^n{ x f 



(2-2) 
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where i is the length of the system, v is the Fermi velocity, a x and a y are the x 
and y components of Pauli matrices, A is the dimensionless coupling constant, loq is 
the phonon frequency and II(x) corresponds to the momentum conjugate to A(x). 
Nonvanishing equal-time commutation relations among those operators are 

{d a (x),d a >(yy} = {e a {x),e a >{y) Jf } = S a yS(x - y) , (2-3) 
[A(x), n{y)] = ih 2 Tr\vLu 2 5(x - y) , (2-4) 

where {^4, B} = AB + BA and [A, B] = AB - BA. As the system is finite, elec- 
tron waves are reflected back at the edges, and the following boundary condition is 
imposed: 

d a (0) = , e a (£) = . (2-5) 

The fields d a and e a correspond to electrons at even and odd sites, respectively, 
of a finite SSH lattice (cf. Appendix A). Note that, instead of d a and e a , the 
original worlPS uses the right- and left-moving electron fields ipR a (x) and ipLa(x), 
respectively, 

V'ftrO) = -^={d a (x) - ie a (x)} , ipLa{x) = —7={e a (x) - id^^x)} . (2-6) 
V2 V2 

The reservoirs are described by 

H B = Y dklhwkLal^a + hw kR bl a b ka } , (2-7) 

where and stand for the annihilation operators of electrons with wave num- 
ber k and spin a in the left and right reservoirs, respectively, and hjj kv {y = 
L, R) are their energies measured from the zero-bias chemical potentials at abso- 
lute zero temperature. The nonvanishing anticommutation relations among them 
are {a^, a\,i a i} = {b^a, b\,, a ,} = 5 aa '5{k — k'). The chain-reservoir interaction is 
assumed to be 

V = J2jdk j^k4(0K CT + hw k di(£)b ka + (/i.e.)} , (2-8) 

where and stand for the coupling matrix elements. 
The Hamiltonian of the whole system is given by 

H = H S + V + H B . (2-9) 

As will be discussed in Appendix A, an open TLM chain described by H corresponds 
to an open SSH lattice that couples with the reservoirs through the end sites and 
the number of whose sites is a multiple of four. 

From (I2-9|) . the lattice distortion A is found to obey the following equation of 
motion: 

= \[A(x,t),H]=n(x,t) 
at in 

^ = = -^o t) + nhvX £ n(x, t)a m 9 a (x, t)j . (2-10) 

Eqs. (|2-9p and (|2-10p are our starting points. 
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§3. NESS Mean-Field Approximation 

In this section, we describe a procedure for evaluating the NESS averages of the 
electron variables for the TLM chain and derive the self-consistent equation for the 
NESS average A(x) = (A(x)) oc of the lattice distortion, which serves as the order 
parameter of the Peierls transition. 

The nonequilibrium steady state under the mean-field approximation is char- 
acterized by (jl-ip where and (5^ are the incoming fields of a^ a and b^ a with 
respect to the mean-field Hamiltonian: 

H MF = Hf F + V + H B (3d) 



„ Jo 



-ihva v -^- + A(x)a x 
y dx 



V ff (x) . (3-2) 



Namely, they are defined as the solution of 

i[a kCT , H U f] = u kL a Vc , jHuvt/H a ^ e -iH Mv t/h e iu kL t ^ ^ (t ^ _ Qo) (3 . 3) 

\\Pk*, H MF ] = uj kR f3 ka , e iH^t/% Kae -iH UY t/n e iu, hR t ^ ^ {t ^ _ oo) _ (3 . 4) 
a 

Since the mean-field Hamiltonian -£Tmf is bilinear with respect to the electron cre- 
ation/annihilation operators, the incoming fields are linear combinations of a^a, b^, 
and ^ CT (x). As shown in Appendix [Bj the incoming fields are fully determined by 
(|3-3p and (|3-4p . Conversely, the original operators can be represented by the incom- 
ing fields. For example, we have 

mr\ Ia\\ H x ^kL) . h(x;u} hR ) \ 

where A-(u>), h(x;uj) and h(x;uj) are auxiliary functions given by 

A_ H = 1 - £_ (co)g^ (0, 0; to) - t?_ {u)g ++ (£, £; u) 

+ £_(u;)ri-(u){g ++ (e,£;u;)g—(0,0;u) - g+^(£,0;u)g^ + (0,£;u)} , (3-6) 

h(x; u) = G(x, 0; u) f Jj {1 - g ++ (£, £; u) V+ (u)} 

+G(x,£;u)r\g + -(£,0;u)r)+(u) , (3-7) 
h(x; u) = G(x, 0; u) Q (0, £; (u) 

+G(x,£;u)r\{l-g^(0,0-,u)C+(co)} , (3-8) 
= / dk> -, r? ± M = / dk' |Wk ''' - . (3-9) 
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And the Green function G for the finite TLM chain and, equivalently, its components 
g aa i (a, a' = ±) are defined as a solution of 

G(x,y;u)=( 9++( ; X > y > U l 9 + -<f,V-,»)\ (3 . 10) 
^ 7 \g-+(x,y;u) g^(x,y;uj)J K ' 



G(x, y : uj) = HujG(x, y; uj) — hl5(x — y) , (3-11) 



0_ 

Ox 

5++(0,y;w) = g + ^(0,y;uj) = g_ + {t,y;u) = g — (£,y;u) = . (3-12) 



-ihvfjy— — h A(x)a x 



For the derivation of incoming-field operators, see Appendix[Bj Then, the mean-field 
NESS is given as a state satisfying Wick's theorem with respect to and (3^ a with 
the two-point functions (jl-ip . Note that, although the Green function (and, thus, 

h and h) diverg function of cj at eigenvalues of the differential operator in 

the left-hand side of (|3-lip . the integrand of (|3-5p remains finite even when uj^l or 
uj^r is equal to one of the eigenvalues. The NESS average of any electron variable 
consisting of & a can be calculated from (jl-ip and (|3-5|) . 

The self-consistent equation for the order parameter A(x) is derived from the 
equation of motion (|2T0p for the lattice distortion. Because of the time-independence 
of A(x), ([2d0]) leads to 

= ^^^=A(x)+7rhv\T(^ a (x)aMx)) oo . (3-13) 

u <j 

Using (jl-ip and ()3-5j) . / dk\v]f\ 2 F{ujki) = / dujF(uj) Im£_(iij)/7r and a similar 

J J — oo 

formula for n~(uj), the self-consistent equation (|3-13|) reads as 



2 -rJ 

^ J-oo 



t c ( Mx;uj)^a x h(x;uj) 
1 Im? -( w ) \A_(uj)\ 2 M > 

h(x;uj)^a x h(x;oj) \ A{x) 

+Imr? " M \A_(uj)\* = "^A ' (344) 



where we use the convention Im£_(w) = (Imi]-{uj) = 0) for uj outside the range of 
^>kL i^kR)- Equation (|3-14p is the self-consistent equation for the order parameter 
A(x). 



§4. Spatially Uniform Phase 



4.1. Self- consistent equation 

Hereafter, we consider the cases where the chain is half-filled; namely, the zero- 
bias chemical potentials are located at the band centre of the TLM chain. Then, in 
order to prevent an increase in electrostatic energy, the chemical potentials of the 
reservoirs should be chosen such that hl = —[J.r = —eV/2 where V is the bias voltage 
and e is the elementary charge (for details, see Appendix [C]). It is well knowrP^J 1 
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that the energy cutoff kw c (Huj c ^> T, e\V\) is necessary for the TLM model, and the 
integration interval of (|3T4p should be replaced with (—u> c ,u> c ). 

The local state of the TLM chain within a certain region of the boundary is 
affected by the existence of the reservoirs. However, if the chain-reservoir interaction 
is not too strong and if the size of this region is much smaller than the length of 
the chain, boundary effects might be neglected. From this observation, we study the 
uniform phase where the order parameter is independent of the coordinate: A{x) = 
A. In this subsection, we derive the self-consistent equation for A when the TLM 
chain is sufficiently long. 

When the order parameter is spatially uniform, the Green function defined by 
(|3Tip can be easily obtained for any real number u and we have (for its complete 
expression, see Appendix [D]) 

hjj sin kx ... 

g ++ (x,e;u)=g—{£-x,0;u) = — — , 4-1 

vL)(uj) 

Hkv cos kx + A sin kx ... 
g- + (x,e;u)=g + -(t-x,0;u) = —- , 4-2 

VJJ{U!) 



where k = y (ftuj) 2 — A 2 /(hv) and D(uj) = Hkv cos k£ + AsitykI. Then, the second 
term of the left-hand side of (|3T4p reads as 

h(x; Lu)^a x h(x; lo) 



f j t / ^ h(x;uya x h{x;u}) 
/ dwlm.r)-{u) MM 

J\u,\<u c l^-Hr 

f hn.T}-(w) f R {hw) r 2 2 2 . 

du - - < FfKvuj(v - \£-{u)) sm2ftx 

M<Wc \v*D(lj)A-(uj)\ 2 I 



- {hwv 2 A + 2v(frw) 2 Re£-(uj) + hujA\^(uj)\ 2 ) cos 2kx 

+ fe 2 Z( W + ^ReC-M + ^|^H| 2 )}sgn(|^| -\A\) , (4-3) 

where the denominator is a function of u and k£: \v 2 D(uj)A_(uj)\ 2 = \A_(lu, k£)\ 2 
with 

A-(u, 0) = fiKv(y 2 — £-(u)rj-(uj)} cos 9 

+ (v 2 A + vhu;{^(u;) +r]^(u)} +A^(co)rj-(uj)) sin9 . (4-4) 



Now we show that (|4-3j) is simplified when the TLM chain is sufficiently long. It is 
easy to see that if | fouj \ < \A\, k is purely imaginary, \D(u))\ 2 h(x; u>)^o x h{x; u) ~ e 2 ' K ' x 
and \D(oj)A_(uj)\ 2 ~ e 2 ' K K Then, the integrand of (|4-3p is on the order of 

h(x;^a x h(x;uj) - 2 \k\(£-x) (a ^ 

\A_(u;)\ 2 ~ B ' 1 J 

which is negligible for large i unless x is near the chain ends. Hence, the contribu- 
tion to (|4-3p from the interval < to < \A\/h is negligible. This implies that \A\ 
corresponds to the electronic energy gap as in the equilibrium case. 
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On the other hand, if \Huj\ > \A\, k is real. Then, with the aid of the Fourier 
expansion of 1/\A-((J, 6)\ 2 , 



where l/£ n (u) is the Fourier coefficient, one can show that the first and second terms 
of the right-hand side of (|4- 3|) are negligible, and that the denominator in the third 
term is replaced with Co( w ) provided I is large and x is far from the chain ends. For 
example, the first term of (|4-3p becomes 



lmr)-(uj) f R (tuj) 2 2 , t / n,2n ■ o 
i o t-,/ \ < 7 \ i o kvluIv — \€-Imj)\ )Sm2KX 
|2|M<M< Wc |^ 2 £(u;)/L»| 2 

where H(u) = lm.rj_(uj) fn(fkj)h 2 kvuj(v 2 — |£_(w)| 2 ). When x is far from the chain 
ends, n^±x = 0(£) and (|4-6p is negligible for large ^ thanks to the Riemann-Lebesgue 
lemma. In short, (|4-3p is found to be 



u)^a x h(x; uj) 

hv 2 A [ ^5^^/ R (M(^ + ^Ree-M + 4ie-MI 2 ) ■ (4-7) 



Jn<co c \ A -M\ 



By a similar argument, when the TLM chain is sufficiently long, the self-consistent 
equation (|3-14|) gives the desired equation for A 

A = or (4-8) 
^ = S(A,V,T L ,T R ) , (4-9) 



where 

S(A,V,T L ,T R ) f Jlm^H / 2A u 2 \ f (h v 

— — = J ^\^7r(" + m7 Rer? - M+ ^ |r? - MI 



J #<l^l<^c 



The function Co( w ) i n the denominators is easily calculated as 



Im£_(cj) <j w + _Re? ? _(a;) + _|^_(a,)|* J. 



+Imr?_(u;) <J w + ^Re£_(u;) + 

ra; ?r 



(4-12) 
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4.2. Electric current 

As easily seen, the electric current at x in the TLM chain is given by 

J(x) = -ev ¥l(x)<T v ¥ (x) , (4-13) 

(7 

and its NESS average by 

-J l7( « ^v 2 p , 2 Im?74^)Im^4^ / \ 

^ = (^))oo = -— J^M*™) \ v 2 D[uj)A _ {uW (/l(M-/MMJ • 

As in the case of the self-consistent equation, for large the averaged current reduces 
to 

J = - 4 -^ [ Mn ^-^-^ - ■ (4-i4) 

^■<\U>\<U1 C 

4.3. Stability 

Since no general thermodynamic criterion is available for discussing the sta- 
bilities of NESS, we study the phase stability based on the linear stability of the 
adiabatic evolution equation for a spatially uniform order parameter. 

Within the adiabatic approximation, the force on the order parameter from 
the electrons is given by -Khv\Y^ ICJ {^ , t{x)o- x ^ a {x)) 00 , where the order parameter is 
replaced with its instantaneous value A(t). Then, we have 



d 2 A(t) 
dt 2 



-co 



{l + \S(A(t),V,T L ,T R )}A(t) , (4-15) 



where the function S(A,V,Tl,T r ) is defined by (|4-10p . 

First, we consider the stability of the normal phase where A = 0. As the 
linearized equation for A(t) is given by 



8 2 A{t) 



§{l + \S(0,V,T L ,T R )}A{t) 



dt 2 

the phase is stable when 

Xn = 1 + XS(0,V,T l ,T r ) >0 . (4-16) 

Next, we investigate the stability of the phase with nonvanishing A. Then, two 
cases should be distinguished: the constant-bias-voltage and constant-current cases. 
In the former, the linearized equation for 5A(t) = A(t) — A is 



and the phase is stable when 

■„(A) = XAl 

dA 



Xv(A) = \A[^=) >0. (4-17) 



v 
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In the latter case, the linearized equation is 



d 2 A(t) J,,/ dS \ ~. , 

-5^ = "^^) J A ®> 



dt 2 

and the phase is stable when 

»o> -*(!)_ 

=°{(a),-(s) 3 (i),./(s)j.- - 

As will be shown in Appendix [Ej the phase with a nontrivial order parameter at 
constant current is more stable than that at constant bias voltage. 

§5. Nonequilibrium Phase Transitions 

5.1. Basic formula 

In this section, we study the nonequilibrium phase transitions when the TLM 
chain couples symmetrically with two identical reservoirs at temperature T: |t>k| = 
kk|, u kL = uj kR , and T L = T R = T. Since r}-{u) = £-{oj) and 



sgn(w) 



for \hui\ > \A\ where sgn((j) the sign of u, we have 

S(A,V,T L ,T R ) =S(A,V,T)= [ — sgnH{/ L (M + /fl(M} • 

T L =T R =T J KV 

^■<\ttl\<W c 

Thus, the self-consistent equation (j4 9j) becomes 

1 arAVT\ ^ de 2sinh(e/r) 

- = -S(A,V,T)= — — . (5-1) 

/ e 2_2\ 2 cosh(|y) + cosh(e/T) 



Since the current is carried by electrons having energies near the Fermi energies, we 
further approximate £-(u;) = rj-(u) = iwo and we have 

J = — rr~r^ 2\ / ^t-t(7l(M - MM) 

^(■U 2 + Wq) _ J \UJ\\ J 



2Go ^ <? - A' sinh(f) 
e e cosh(|^) + cosh(^) 
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where the normal-state conductance Gq is given by 



Go 



2vwq 



ttH v 2 + Wq 



(5-3) 



Indeed, if A = and terms on the order of exp{— {2fno c — e|V|)/(2T)} are neglected, 
one can easily evaluate the integral in (|5-2p . yielding 



J = GqV , 

irrespective of the temperature T. In the rest of this section, the phase transition 
and nonlinear conduction will be discussed based on (15-11) and (15-21). 



5.2. Phases at absolute zero temperature 

In this subsection, phases at absolute zero temperature are investigated. Thanks 
to the formula 



sinh(y/r) 
T^o cosh(y/T) + cosh(x/T) 



sgn(y)6(\y\ - \x\) 



with the step function 9, the function S(A, V,T) and the current at T = are given 
by 



S(A,V,0) 



hjJc 



de 



Ml 



eV 



(5-4) 



J = — sgn(y) / de 
e Jm 



i-A z 



eV 



The stability index of the normal phase A = is, then, 

I VI 

XJV = l + AS(0,y,0) = 2Alni^ , 

no 



(5-5) 



(5-6) 



where V\q is defined by V\q = ^|^exp(— Jv-). Hence, the normal phase is stable if 
|V| > Vio and unstable if |V| < V\q. As mentioned in the previous subsection, the 
average current is given by 

7 = GqV 
For |Z| > |eV|/2, ([53]) reduces to 

"JjjJJc 



S(A,V,0) 



de 



Ml 



e 2 -A 2 



-2 cosh 



huj c 



(5-7) 
(5-8) 



and (|5-ip has a nontrivial solution, 

01 



cosh^ 



(5-9) 
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irrespective of the bias voltage \V\ < V20 = 2Aq/c. Its stability index at constant 
bias voltage is always positive: 

/dS\ 2\hw c 

XV = A^(^) v = . > . (5-10) 



In short, the ordered phase exists for | V| < V20 and is stable at constant bias voltage. 
According to the inequality xi > XV shown in Appendix lEj this phase is also stable 
at constant current. Since the Fermi energies of the two reservoirs fall into the energy 
gap, the phase is insulating: 

7 = 0. (5-11) 

Because of V20/V10 = 2/(1 + e" 1 ^) > 1, there might be a first-order phase 
transition between the normal and insulating (\A\ = Aq) phases for V10 < \V\ < V20. 
This suggests the existence of another solution \A\ of (|5-ip satisfying < \A\ < Aq. 
Indeed, for |Z| < |eV|/2, 



S(A,V,0) = -2 ^ ; d€ = -2{ cosh- 1 cosh" 1 S) , (5-12) 



rftu) c 

and thus (|5-ip has a nontrivial solution for V10 < |V| < V20: 




This phase is unstable at constant bias voltage since the stability index is negative 
for \A\ < \eV\/2 < frw c : 

,-rfdS\ 2\hw c \\eV\ n , ri/1 , 

Xv = AZ\ ^ = — _ 1 1 == < . (5-14) 

KdAJv J{hw c y-\A\* J{eVf/±-\A\> 



In contrast, at constant current, this phase is stable because the stability index is 
positive: xi > (for a proof, see Appendix [E]). In this case, the average current is 
given by 



2G n f^ 2 , \/e 2 -A 2 



- 2G nn [W, 
J = sgn(F) / de 

e J\A\ 



2G> 



1^1 

e ° |Z|sgn(V){ y (g) 2 - 1 - tan" 1 {^f - l} = J 2 (X,V) . (5-15) 

These results are summarized in Figs. QJH At constant bias voltage (cf. Fig. Q]), 
the insulating phase is stable up to the first threshold voltage eVio, which is almost 
half of the zero-bias gap 2Aq for small A. Beyond the second threshold voltage eV2o 
which is equal to the zero-bias gap 2Aq, only the normal phase is stable. Between the 
two threshold voltages, i.e. Viq < \V\ < V20, both the insulating and normal phases 
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are stable, and a first-order phase transition between them is possible. Moreover, 
there exists an unstable phase separating the two stable ones, as shown by the 
dashed curve (i.e., the curve satisfying \A\ < e\V\/2) in Fig. [TJ The current-voltage 
characteristics are shown in Fig. [2j Note that the first-order transition corresponds 
to a sudden change in the current. The three regions discussed above are summarized 
in Table E 

On the other hand, at constant current (see Fig. [3|), the unstable phase men- 
tioned above is stabilized, and when one increases the current, a second-order phase 
transition to the normal phase occurs at the critical current \J\ = J c0 = Go^io- Near 
the critical current, the order parameter changes linearly with respect to the current: 

— 2A J c0 - | J\ 

The corresponding voltage-current characteristics are shown in Fig. a region < 
J < J c o with negative differential conductivity appears. 

The existence of negative differential conductivity can be understood as follows. 
When no bias voltage is applied, the system is in an insulating phase with an energy 
gap 2Aq, which is robust against low bias voltages. When the current starts to flow, 




Fig. 1. Bias- voltage dependence of the order 
parameter at absolute zero temperature. 
In the constant-bias-voltage case, the solid 
lines correspond to the stable phases and 
the dashed curve to the unstable phase. All 
phases are stable in the constant-current 
case. 



Table I. Three regions at absolute zero temperature (constant-bias-voltage case) 





V 


\A\ 


J 


region A 


< V < Vio 


A)(A) 





region B 


v 10 <v< v 20 


0, A (X), Z\i(A,V) 


G V, 0, J 2 (X,V) 


region C 


v 20 <v 





G V 



Fig. 2. Current versus bias voltage at abso- 
lute zero temperature. The solid lines and 
dashed curve correspond to those of Fig. [1] 
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< 1 

l<3 



0.8 




0.5 1 1.5 2 _ 2.5 

J/J cO 



Fig. 3. Current dependence of the order pa- 
rameter at absolute zero temperature. In 
the constant-current case, all phases are 
stable. Note that the point (J, \A\) = 
(0, Ao) corresponds to the insulating 
phases for Vio < |V| < Vao- 




0.5 



0.5 1 1.5 2 _ 2.5 

J/J cO 



Fig. 4. Bias voltage versus current at absolute 
zero temperature (constant current). The 
curve is essentially the same as Fig. [2] Only 
stabilities are different. 



electrons with energy larger than the gap should exist, and the bias voltage should 
be on the order of the zero-bias gap: \V\ ~ 2Ao/e. At the same time, the gap is 
reduced by the existence of the current. As the current increases, the gap 2|Z\| is 
reduced further, and the corresponding bias voltage |V| ~ 2\A\/e becomes smaller. 
Thus, negative differential conductivity does appear. When the current reaches the 
critical current, the gap disappears and the phase becomes normal. As the normal- 
phase conductivity is positive, negative differential conductivity appears only up to 
the critical current: |J| < J c q. 

5.3. Phases at finite temperature at constant bias voltage 

In this subsection, we investigate the finite-temperature phases at constant bias 
voltage. 

(A) Phase Diagram 

Let us begin with the investigation of the phase diagram. As in the zero- 
temperature case, three regions exist in the l^T-plane; one with a unique stable 
ordered phase (region A), one where ordered and normal phases are stable (region 
B) and one with the normal phase (region C). The three regions are depicted in 
Fig. 5. The boundary curve between region A and the others is implicitly given by 

1 l„,„„„s 1 /"^ c/T de sinhe 



2A 2A 2 y 1 2A J 



cosh(|K) + cosh e 



~2T' 

where xn is the stability index of the normal phase. When terms on the order of 
exp{ — (2hui c — e|V|)/(2T)} are neglected, it reduces to 



eVw f°° , , cosh ( f^) cosh e + 1 jeV 

= / de loge % = 61 — 

2T J & (cosh(f^) +coshe) 2 "v2T 
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where the function (ft is defined by the integral of the middle term. When V 
(fFTTj) leads to 



log 



2T 



de 



loge 



1 + cosh e 



log 



2e^ 



7T 



or 



-cO 



-eVi 



7T 



10 



(5-18) 



where 7 is the Euler constant and T c q corresponds to the transition temperature at 
zero bias voltage. Then, (I5-17|) reads 



log 



/ 7T T c0 
\2e-y T 



A2eTVio T 



(5-19) 



This implies that the boundary between region A and the others is independent 
of the coupling constant once it is plotted in terms of V/V10 and T/T c q (cf. the 
solid curve in Fig. [5]). We note that the boundary curve (|5-19p can be expressed as 
|V| = V\{T) in terms of a single-valued function V\{T) of T, which will be referred 
to as the first threshold voltage. 

The phase boundary curve (|5-19p indicates the bias-induced decrease of the 
transition temperature for low \ V\ and the temperature-induced increase of the first 
threshold voltage for low T. Indeed, for e\V\ <C T, with the aid of the formula 
(ft'(0) = and 



0"(O) 



00 ^ loge(coshe - 2) 



(1 + cosh e)' 



7C(3) 
2vr 2 



(|5-19p reduces to 
T = T C (V) 



! c0 exp 



7C(3) (VT c0 y 
IGe 2 "/ \V 10 TJ 



* Tcol 1 



7C(3) / V n 
16e 2 ^ \V W 



) 2 } , (5-20) 



where C( n ) is the Riemann zeta function. Thus, the transition temperature T C (V) 
decreases with the bias voltage V. On the other hand, when T -C e\V\, one may 
apply the standard technique of evaluating the low-temperature properties of the 
free fermion gas, and we have 



delog (-^=) — 
B \2TJ de 



1 



+ 



1 



(e-eV)/(2T) + 1 e (e+eV)/(2T) + j 



log 



e\V\ 
2T 



2tt 2 /T\2 



log 



e\V\ 2e 2 "i (TV 10 \ 2 



2T 



3 \T c0 VJ 



(5-21) 



up to T 2 . Thus, the first threshold voltage V\{T) near absolute zero temperature is 
given by 



|V| =V 1 (T) -Fioexp 



2e 2 ^/TV 10 y 

3 vr c0 Vi/ 



10 



2e 2 ^ / T \ 

1 + — 



} . (5-22) 



which shows that the first threshold voltage increases as temperature increases. Care- 
ful asymptotic analysis indicates that T/(e|V|) ~ 0.05 is the upper bound where the 
estimation (15-22p is valid. 

On the other hand, the boundary curve \V\ = V^(T) between regions B and C 
(the dashed curve in Fig. [5]) is derived by solving the self-consistent equation (|5-ip 
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for nonvanishing order parameters. It starts from the point (|V|,T) = (V2o,0) and 
terminates at a point P on the boundary curve | V| = V\{T) between region A and the 
others (see Fig. 5). The behaviour near the terminating point P can be investigated 
on the basis of the Ginzburg-Landau expansion of the self-consistent equation (|5T|) . 

— 4 

Up to A , the self-consistent equation becomes 



2 \t) 8 \t) 2a 



(5-23) 



where xn = 1 + AS(0, V, T) is the stability index of the normal phase, and the 
coefficients Ki and K± depend only on the ratio V/T (their concrete expressions 
are given in Appendix |F|). As V/T 
increases, K2 changes sign from mi- 
nus to plus at V/T ~ 2.1865 x V 10 /T c0 
and K4 is positive there. Hence, when 
K2 < 0, the quartic polynomial in the 
left-hand side of f)5 • 23|) has one max- 
imum value at \A\ = 0. In this 
case, when xn > 0, f)5 • 23|) has no so- 
lution (region C), and when xn < 0, 
(|5-23p has one nonvanishing solution 
in \A\ (region A). On the other hand, 
when K2 > 0, the quartic polyno- 
mial in the left-hand side of (I5-23|) 
has a local minimum at \A\ = 0, 
and a local maximum K\/(2K^) at 
|Z| = Ty / 2K 2 /K i . Then, with re- 
spect to |Z\|, (|5-23p has no solution 
when K$/(2K 4 ) < Xtv/(2A) (region 
C), two nonvanishing solutions when 
Kl/{2K A ) > xn/(2X) > (region B) 
and one nonvanishing solution when 
Xtv/(2A) < (region A). Therefore, 
the boundary curve \V\ = V^(T) near 
point P is given by K%/(2K±) = 




Fig. 



5. Phase diagram at constant bias voltage. 
A unique stable ordered phase occurs in re- 
gion A, one stable ordered phase and the nor- 
mal phase in region B and only the normal 
phase in region C. The solid curve represents 
the second-order phase transition, and the 
dashed curves represent the first-order phase 
transition. For comparison, the phase transi- 
tion curve for (a) A = 4.8 is shown together 
with that for (b) A = 3.0. 



Xjv/(2A). Then, since the curve \ V\ = 

V\(T) corresponds to xn = 0, the simultaneous solution of (|5- 19j) and K2 
the terminating point P: T = T* ~ 0.5571 x T c0 and \V\ ~ 1.2181 x Via. 



-- is 
Note 

that, in contrast to the first-threshold-voltage curve |V|/Vio = Vi(T) /V\o, the curve 
|V|/Vio = V2(T)/Vio depends not only on T/T c q but also on A (cf. Fig. [5]). 



(B) Order Parameter 

Both the bias- voltage dependence of the order parameter and the current- voltage 
characteristics change continuously from those at zero temperature T = as shown 
in Figs. [6]l9] for 1/A = 4.8. When temperature is less than T* corresponding to 
the point P of Fig. 5, a voltage range Vi(T) < \V\ < V^(T) exists where the order 
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0.2 



1 1 i i i i I 

0.5 1 1.5 2 2.5 

V/Vw 



Fig. 6. The order parameter versus bias volt- 
age for (a) T = 0.01exp(-l/2A)|eVio| (= 
0.194T c0 ) and (b) T = 0.02 exp(l/2A)|eVio | 
(= 0.389T c o). The solid lines correspond to 
the stable phase, and the dashed lines cor- 
respond to the unstable phase. 

a 1.2 , , , , . 

< 




V/Vw 




0.5 1 1.5 2 



Fig. 7. Current versus bias voltage for (a) T = 
0.01 exp(l/2A)|e Vio| (= 0.194T c0 ) and (b) 
T = 0.02exp(l/2A)|eVio| (= 0.389T c0 ). 
The solid lines correspond to the stable 
phase, and the dashed lines correspond to 
the unstable phase. 




0.5 1 1.5 2 



V/Vw 



Fig. 8. The order parameter versus bias volt- 
age for (a) T = 0.03exp(l/2A)|eVio| (= 
0.583T c0 ) and (b) T = 0.04exp(l/2A)|eVio| 
(= 0.778T' c o). All phases shown in this fig- 
ure are stable. 



Fig. 9. Current versus bias voltage for (a) T = 
0.03exp(l/2A)|eVi |(= 0.583T c0 ) and (b) 
T = 0.04exp(l/2A)|eVio|(= 0.778T c0 ). All 
phases shown in this figure are stable. 



parameter is a triple- valued function of the bias voltage |V| (see Fig. [6]) with an 
unstable middle branch and, hence, where the first-order phase transition is possible. 
The corresponding current-voltage characteristics are S-shaped, as shown in Fig. [71 
Note that the small current observed at low bias voltage is due to thermally activated 
carriers. On the other hand, when the temperature is higher than T*, the unstable 
branch disappears and the order parameter becomes a single-valued function of the 
bias voltage (see Fig. [5}. In this case, the current is a monotonically increasing 
function of the bias voltage (see Fig. [9|) . 
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< 1 



I<1 




0.2 0.4 0.6 0.8 1 

T/Tco 

Fig. 10. The order parameter versus tempera- 
ture for V = 0.551Vio. The solid lines cor- 
respond to the stable phase and the dashed 
lines correspond to the unstable phase. 
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< 
I<1 




0.2 0.4 0.6 0.8 1 



T/Tco 

Fig. 11. The order parameter versus tempera- 
ture for V = 1.047Vio- The solid lines cor- 
respond to the stable phase and the dashed 
lines correspond to the unstable phase. 




Fig. 12. The order parameter versus tempera- 
ture for V — 1.213Vio. The solid lines cor- 
respond to the stable phase and the dashed 
lines correspond to the unstable phase. 



Fig. 13. The order parameter versus tempera- 
ture for V = 1.224Vio- The solid lines cor- 
respond to the stable phase and the dashed 
lines correspond to the unstable phase. 



The temperature dependence of the order parameter at constant bias voltage 
is shown in Figs. [T0II151 At low bias voltage, the temperature dependence of the 
order parameter is similar to that in the absence of the bias (see Fig. [TUj) . As 
the bias voltage increases, a lower-temperature branch corresponding to the un- 
stable phase appears (see Fig. ITT]) . In the temperature range where the unsta- 
ble phase appears, the normal phase is stable in the sense of xn > 0, and the 
first-order phase transition between the stable ordered and normal phases is pos- 
sible. As the bias voltage increases further, the unstable branch approaches the 
stable branch and the two branches join (Figs. fT2lll4|) . Note that an unstable por- 
tion appears in the outer curve of Fig. [13] This is because the temperature T* is 
higher than the temperature where = 0. For higher bias voltage, the sta- 

ble ordered phase always co-exists with the normal phase, and the region where 
they co-exist shrinks with increasing bias voltage (see Figs. [TJ] and [T5]) . We re- 
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Fig. 14. The order parameter as a function of 
temperature for V = 1.235Vio. The solid 
lines and dashed lines in this figure have 
the same meaning as those in Figs. [6] and 

Fig.m 



Fig. 15. The order parameter as a function of 
temperature for V = 1.433Vio. The solid 
lines and dashed lines in this figure have 
the same meaning as those in Figs. [6] and 

Fig.m 



mark that both the phase diagram (Fig. [5]) and the temperature dependence of 
the order parameter (Figs. Q3] and [T5]) are similar to those for the nonequilib- 
rium superconducting phase induced by excess quasiparticles, which was studied 
by Scalapino et alW^^b This will be discussed in more detail in the last section. 



5.4. Phases at finite temperature at 
constant current 

(A) Phase Diagram 

As in the zero-temperature case, 
at constant current, all the nontriv- 
ial solutions are stable in the sense 
of xi > 0] an d the properties are 
drastically changed from those at con- 
stant bias voltage. First, only the or- 
dered and normal phases exist and the 
phase transition between the two is 
always of second order. The corre- 
sponding boundary curve is given by 
{FIJI: 



log 



ttT c qJ 



H 



0.6 



0.4 



0.2 





Normal 


Ordered 





0.2 



0.4 



0.6 



1.2_ 1.4 
J/J cO 



Fig. 16. Phase diagram at constant current. 
Only the second-order phase transition occurs 
at constant current. 



\2e~<TJ T \2e~tJ c0 T J 

(5-24) 

where J c n is the zero-temperature 

threshold current defined just above (|5-16p (Fig. I16[) . From (|5-24p . the critical tem- 
perature T C (J) for small current and the threshold current J C {T) at low temperature 
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0.2 0.4 0.6 0.8 1 1.2_ 1.4 



J/J cO 



Fig. 17. The order parameter versus current 
for (a) T = 0.01exp(l/2A)|eVio| (= 
0.194r c0 ) and (b) T = 0.02 exp(l/2A)|eVio | 
(= 0.389T c o). All phases shown in this fig- 
ure are stable under constant current. 

are found to be 




0.5 



0.5 1 1.5 _ 2 

J/J cO 



Fig. 18. Current versus bias voltage for (a) 
T = 0.01exp(l/2A)|eVio| (= 0.194T c0 ) and 
(b) T = 0.02exp(l/2A)|eVi | (= 0.389T c0 ). 
All phases shown in this figure are stable 
under constant current. 



and 

\J\ = J C (T) ~ J c0 {l + -3"(j^) } • (5-26) 

The difference between the phase diagram at constant bias voltage and that 
at constant current can be understood as follows. Since a larger order parameter 
implies a smaller current, the nontrivial phase with larger order parameter in region 
B of Fig. [5] corresponds to the phase with smaller current. As a result, the phase 
diagram in Fig. [16] does not have a region where more than one phase is stable. 

(B) Order Parameter 

As before, when temperature increases, both the current dependence of the 
order parameter and the voltage-current characteristics change continuously from 
those at zero temperature. At any temperature, the order parameter is a mono- 
tonically decreasing function of the current (cf. Figs. [T71 and I19p . For |J| ~ 0, the 
decrease in the order parameter is found to be proportional to the squared current, 
but the quadratic region is not visible at lower temperature T < 0.6 x T c q ~ T* , 
where the decrease in the order parameter is approximately proportional to the cur- 
rent. The corresponding voltage-current characteristics behave as shown in Figs. [18] 
and [20l We remark that negative differential conductivity appears only at tem- 
peratures lower than T*, and otherwise the differential conductivity is positive. 
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J/J cO 

Fig. 19. The order parameter versus current 
for (a) T = 0.03exp(l/2A)|eVio| (= 
0.583T c0 ) and (b) T = 0.04exp(l/2A)|eVio| 
(= 0.778T c o). All phases shown in this fig- 
ure are stable under constant current. 



The temperature dependence of 
the order parameter at constant cur- 
rent is shown in Fig. [2T1 The or- 
der parameter is reduced by the pres- 
ence of the current. Even for a 
smaller current (see the cases J / J c q = 
0.3 and 0.6 in Fig. [21]), the lower- 
temperature part is more suppressed 
than the higher-temperature part. At 
the threshold current J = J c q, the 
order parameter at zero temperature 
vanishes (see the case J / J c q = 0.95 
in Fig. I2T1) . When the current ex- 
ceeds the threshold, re-entrance to the 
normal phase appears at low temper- 
ature. As the current increases, the 
temperature range with nonvanishing 
order parameter shrinks and eventu- 
ally vanishes (see the case J/J c o = 
1.10 in Fig. EH). 

§6. Summary 
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J/J cO 



Fig. 20. Current versus bias voltage for (a) 
T = 0.03 exp(l/2A)|e Vio| (= 0.583T cO ) and 
(b) T = 0.04exp(l/2A)|eVio| (= 0.778T c0 ). 
All phases shown in this figure are stable 
under constant current. 




T/Tco 



Fig. 21. The order parameter as a function of 
temperature for J / J c o = 0, 0.3, 0.6, 0.95, and 
1.10 from top to bottom. To draw this figure, 
we use an approximate equation of (|5-1[) sim- 
ilar to (|5T7p . For computational reasons, the 
numbers of data are limited for J/J c o =0.3, 
0.6, 0.95 and 1.10. 



and Discussions 



We have studied the nonequilibrium Peierls transition in a TLM chain connected 
to two reservoirs at different chemical potentials (their difference corresponds to the 
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bias voltage) by combining a mean- field approximation and formula which 
characterizes a nonequilibrium steady state and is an outcome of the algebraic field- 
theoretical approach to nonequilibrium statistical mechanics. The averaged lattice 
distortion serves as an order parameter, and its self-consistent equation is obtained by 
averaging the equation of motion of the lattice distortion with respect to a nonequi- 
librium steady state. When the bias voltage and temperature are chosen as control 
parameters, three parameter regions are distinguished: region A where a single stable 
ordered phase is possible, region B where stable normal, stable ordered, and unstable 
ordered phases are possible and region C where only the stable normal phase is pos- 
sible (cf. Fig. 5). The transition between regions A and C is of second order, and the 
transition temperature decreases with increasing bias voltage. A first-order phase 
transition between normal and ordered phases may occur in region B. In regions A 
and C, the current is a single-valued function of the bias voltage, and in region B, 
the current-voltage characteristics are S-shaped (i.e. negative differential conduc- 
tivity exists). In contrast, when the current and temperature are chosen as control 
parameters, all the nontrivial solutions of the self-consistent equation become stable 
in the sense of linear stability. The phase transition between the ordered and normal 
phases is always of second order, and re-entrant behaviour is seen for currents larger 
than the threshold value, J c q. Negative differential conductivity appears only when 
the temperature is lower than a certain value T* . We remark that, as in the equilib- 
rium case,® the mean field approach is expected to provide a qualitatively correct 
description of the nonequilibrium phase transition in quasi-lD systems, although the 
ID order may be destroyed by fluctuations. 

As mentioned in the previous section, the phase diagram, particularly the pos- 
sibility of the re-entrant behaviour, and the temperature dependence of the order 
parameter at higher bias voltage are similar to those of the nonequilibrium supercon- 
ducting phase induced by excess quasiparticles, which were studied by Scalapino et 

al m\m 

This can be understood because of the similarity between the self-consistent 



equations in the two cases. Indeed, let E be y e 2 + A , then the self-consistent equa- 
tion (|5-ip is rewritten as 



which reduces to the self-consistent equation of Scalapino et al. [cf. Eq.(6) of Ref. l32p ] 
if the second term is dropped and eV/2 is replaced with the effective chemical po- 
tential fj,* . In view of this similarity, one can interpret the suppression of the charge- 
density-wave order induced by the bias voltage (equivalently by the current) as being 
due to excess electrons coming from the two reservoirs. However, since the two self- 
consistent equations are not exactly the same, the two systems are different in their 
properties at lower bias voltage (namely, in the temperature dependence of the order 
parameter and the possibility of the second-order phase transition) . 

As is well known, systems with density waves may exhibit nonlinear conduction 
due to density-wave sliding.® On the other hand, because of the above observation 
and the similarity between the mean-field approximation for superconductors and 
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that for density waves, the current-induced suppression of order discussed here is 
generally expected for systems with density waves. However, except for cases where 
density waves are strongly pinned, the current-induced suppression of order may not 
be observed because it is related to the amplitude degrees of freedom, while the 
sliding is related to the easily excitable phase degrees of freedom. 

We remark that the mean-field approximation of the TLM chain is equivalent to 
the mean-field approximation of the half-filled charge order in the spinless extended 
Hubbard chain, the Hamiltonian of which is given by 

L-l L-l L 

Hco = ~t £{Ct +1 Cj + C]C j+1 } + rij+mj -U^nj , (6-2) 

3=0 j=0 j=0 

where rij = C^Cj is the number operator of the spinless fermions at site j. Indeed, 
in the continuum limit discussed in Appendix A [cf. (|A-2[) ] . the mean-field equations 
become (|3-2[) and (|3-13p with replacement a x — ► a z and A — > U/(irto), where the 
order parameter is proportional to the charge disproportion between even and odd 
sites. Hence, as in the open TLM chain, the charge order in the open extended 
Hubbard chain is suppressed by current. This observation suggests that current- 
induced suppression could be a possible origin of the nonlinear conduction, which is 
different from phenomena such as sliding density waves strong impurity scattering 
in a Tomonaga-Luttinger liquidp^ dielectric breakdown of Mott insulators®} or the 
Kosterilitz-Thouless transitionPSJ 

It is then interesting to compare the present results with the experiments by 
Terasaki and his co-workers^J'^'^SJ'^B Q n the charge order in the organic con- 
ductors 0-(BEDT-TTF) 2 CsM(SCN) 4 (M =Zn, Co, Coo. 7 Zno. 3 ). The existence of 
hysteresis and spontaneous oscillation at constant bias voltage and negative differ- 
ential conductivity at constant current reported in Ref. |25|) is consistent with the 
present results. In Ref. |27|) . it was shown that the current-induced decrease in the 
order parameter is proportional to the current and, this observation agrees with the 
current dependence of the order parameter at temperatures where nonlinear con- 
duction is possible (cf. Figs. [T71 and [T9|) . Also, according to the present analysis, 
negative differential conductivity appears only when A ~ U/(irto) is not too large. 
This seems to imply that systems with weaker charge ordering are more favourable 
for negative differential conductivity, and further imply that, because of the fragility 
of their charge order, the organic conductors mentioned above would be such sys- 
tems. Moreover, in these materials, the charge order has no long-range order in the 
phase mode, and thus excitations in the amplitude mode instead of charge-order 
sliding would be responsible for nonlinear conductivity. These observations suggest 
that current-induced suppression of the charge order is the origin of the phenom- 
ena investigated by Terasaki and co-workers.^ More quantitative analysis 
based on a realistic model of the organic conductors will be studied elsewhere. 

Moreover, since the extended Hubbard model (|6-2|) is a typical model of strongly 
correlated systems, the present analysis would provide some insight into the negative 
differential conductivity recently reported in strongly correlated systems.^ JUSMHOJUD 
In particular, as the extended Hubbard chain (|6-2p is equivalent to the XXZ chain 
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via the Jordan- Wigner transformation, the negative differential conductivity in the 
nonequilibrium XXZ model found by Benenti et alW^ should be approximately un- 
derstood in terms of the present results. However, further investigation is necessary 
because in the work of Benenti et al, the system is driven to a nonequilibrium 
steady state by stochastic activation of the boundary spins, not through coupling 
with infinitely extended reservoirs. 
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Appendix A 

Derivation of a continuous model 



In this appendix, we derive the continuous open TLM model from the discrete 
open SSH model. Hamiltonian of our system is consist of the SSH part (H$), two 
reservoirs (Hb), and the interaction between SSH part and reservoirs (V): 

H = H s + V + H B 

L+1 K L M L 

H s = -J2Y1 (wi,nCt +1<T c nCT + (h.c.)) + ^Y1 - y^ 2 + yE^ 

<j n=— 1 n=— 1 n=0 

V = ^2 / dkh (v k Cl a a ka + w k C\ a b ka + (h.c.)) , 
H B = y^l dk (hw kL al a a ka + hv kR b\ a b ku ) , (A-1) 

cr J 

where C na denotes the annihilation operator of an electron at the nth site with 
spin a (cf. C^\ a = 0, Cl+i<j = 0), a ka (b ka ) denotes the annihilation operator of 
the left (right) reservoir with wave number k and spin a, and y n denotes a lattice 
displacement of the nth site, respectively. Su, Schrieffer, and Heeger^J assumed that 
t n +i tn is a linear function of the lattice displacement: 

t n +i,n = to - a(y n+ i - y n ) . 

Takayama et alW^ approximated the dispersion relation of electrons as — 2to cos[(/c± 
&i?)a] w ±vjk, and introduced the left/right moving electron fields ■0Lo-(2na)/'0R fT (2na), 
and A(na) = (— l) n 4ay n with a the lattice constant. Then, by assuming L = —1 
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(mod 4) and considering a to be very small, the discrete Hamiltonian (|A-lj) reads as 



„ Jo v 



d 

-2iatoa y — + a x A(x) 



H 



(Ph) 



K 

8a 2 a 



dxA(x) 2 + dx 



M 



32a 2 a 



A{xf 



H B = y^l dk (hw k Lal a a k a + ^kRb ka ha) ■ 

a J 

Moreover, C-\ a = and Cl+Io- = lead to 

^La{0) + iVRa(O) = 0, i^ La (l) + 4> Ra {l) = , 

where I is defined by £ = a + La. We remark that case of L = +1 (mod 4) leads to 
essentially the same results. To make the boundary condition simpler, we introduce 
d a (x) and e a (x) by 

C 2n -ia = {-l) n V2a~d a (2na) = {-l) n ^[^ La {2na) + i^ Ra {2na)] 

C 2 na = {-l) n V2a~e a (2na) = (-l) n v^frM2na) + ^(2na)] . (A-2) 

By using these fields, the momentum il(x) = Z\(x) conjugate to Z\(x), the Fermi 
velocity v = 2ato/h, dimensionless coupling constant A = 4a 2 a/irhvK , phonon fre- 
quency u>o = y/4K/M, and the matrix elements v k = \[2av k and w k = V2aw k , 
the Hamiltonians (|2-2p . (|2-7p . and (|2-8p . and the boundary condition (|2-5p in §2 are 
obtained. 

Appendix B 

Normal modes 



In this appendix, we derive normal modes of the mean-field Hamiltonian Hmf- 
Let {4>\( x )}\ be a complete orthonormal solution of the eigenvalue problem: 



d .. . 

-lhva v — h ZilXliTa; 

y dx 



(j>x(x) = he x 4>\(x) 



M^^i^X ^(0) = 0, ^W = o 



and we expand the electron field 



d a (x) 
e a (x) 



(B-l) 
(B-2) 

(B-3) 



A.c 



where {/ao- ; = 5a,a"^o-,<t'- In terms of /a ct , the mean-field Hamiltonian reads 
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+ E/ dk {(^( )<i + ^W<t)/^ + M} • 

A,(T 

Since Hup is bilinear with respect to field operators, the incoming fields are linear 
combinations of ako-^ko- and f\ a : 

a-ka = «.k<r + ^2 ^a/ao- + / dk' (m*,av a + n£,ft k / CT ) , (B-4) 
A •* 

By substituting it into [a^, Hyip\/H = uJkL&ka and comparing term by term, one 
has 



Oka = Oka + 2^ "A/Aa + dk [ - — - + 



(0) 



UkL — £A WfcL — £A 

where A£(x) = £ A ^aO)*^ (p = ±; s = 0, and 



(B-5) 
(B-6) 



dk' 



, »7±(*) 



dk' 



(B-7) 



Substituting (|B-6P into the definition of j4^(x), one obtains a linear equation for 
^4 k (0) and A^(£) and its solution is 



A+(£) = v 

«£ + e±(«fcL)^(o)=«j; 



k ^±(^l) ' 
1 - r]±(LU kL )g ++ {l,l : u; kL ) 



(B- 
(B- 



where 



4t(*) = 1 - £t(*)0— (0, 0; z) - rj ± (z)g ++ (£, £; z) 

+ ^±{z) V± (z){ 9++ (£,£;z)g^(0,0;z)-g+^(£,0;z)g- + {0,£;z)} (B-10) 



and g pp '(x,y : z) is the p/Z-component of the Green function: 

'g ++ (x,y;z) g + -(x,y;z)\ _ ^x{x)4>x{y)^ 



G(x,y;z) 



9-+{x,y;z) g — (x,y;z) 



z - e x 



(B-ll) 



As a result of the completeness of the eigenfunctions cpx(x), the Green function 
satisfies 



-ihvcFy— — h A(x)a x 
ox 



G(x, y : z)= hzG(x, y; z) - hl5(x - y) , (B-12) 
9++(0,y;z) =5f+_(0,y;z) =g^ + {£,y;z) =g—(£,y;z) = , (B-13) 
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where 1 stands for the 2x2 unit matrix. Similarly, we have 

a. = + + U ( jvSl^ + -vqrw^j (B . 14) 

h\ = 0a(Q) ^(0)C±Kii) + {< + B+(^( WfcR )} , (B-15) 

Wfci? — €a io kR - e\ 

where B^(x) = ^ A <j)^(x)*h\ and 

Bk(0) = u;k ' (B ' 16) 

< + r ?± (^)i ?k W = ffl£ l-e fc (^-(0,0:o, fcfl ) 

Now we discuss the sign of small imaginary parts in the energy denominators. 
The sign should be chosen so that we have e lHmvt / h a\ L(T e~ lHm¥t / n e %u)kLt — > and 

original 

operators can be expressed in terms of the incoming fields and e.g., 

e iH UF t/h a e ~iH MV t/h e iw kL t _ a ^ 

= v * f dk , |- ^(0)*« k ye-^^^) f | j^,(0)*/W~ tK, *-" fcL)f -( 

k J < CJk'L - ^kL ± ^0 UJk'R - UJkL ±iO > 

which vanishes as t — > — oo only if the lower sign is chosen since lim — = 0. 

t— >-oo x + iO 

Then, the electron field in the TLM chain is given by 

A A 

q^ t( ' ; " tl) ^-h,)- ' l(i; °'' ) } • (B,18) 



where 



h(x;cj) = G(x,0;w)rj{l - <?++(£, 4 

+G(x,£;a;)Qj^ + _(^0;a;)»7 + (w) , (B-19) 
fc(x;o;) = (?(i l 0;«;)rJj- + (0 l <;( 1 )^H 

+G(x,£;o;)Q{l- 5 __(0,0;a;)C+H} . (B-20) 

Appendix C 

Coulomb energy and chemical potentials 

Applying a bias voltage V to the TLM chain corresponds to the prescription 
— fJ-R = —eV and it alone does not determine individual values of fj,i and [Ir. 
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However, as will be explained below, if fi^ + [i R ^ 0, the average number of electrons 
on the TLM chain increases or decreases and the whole system including ions would 
be electrically charged as compared with the equilibrium case. Then, such states 
have large electrostatic energy and are hard to be realized. Thus, one should choose 

HL = -fJ>R = -eV/2. 

The proof is as follows. By a similar argument to the calculation of S in (14- lOj) . 
the electron number density is found to be 



[x) , 

oo 



E(^o^ 

a 

= / -=- i (h(e) + h(-e) + Me) + f R (-e)) , (C-l) 



provided that x is not close to the chain ends. At equilibrium where = ^R = 0, 
the sum of four Fermi distribution functions is equal to two and, thus, irrespective 
to the temperature, 



^ de 2e 



£(*t(s)!M*)) =/ ^- t -^=. (c- 2 ) 



e 2 -A 



Then, as easily seen, one has 



oo * — \ / eq 

u a 

= (1 - e -(^*)/T) de ^)f R {-e) + f L {-e)f R {e)) 

h TThv \je 2 -A^ 

which is nonzero unless hl = —piR. 

In the rest of this appendix, we show that a number of electrons per site for 
the open SSH chain discussed in Appendix A is approximately unity by choosing 
ftw c = 7rfo- Indeed, since 

1 f^" . £ LO r 



E <4(*)*r(^)> = E <4(*M*)> — / d6 



and v = 2ato/h, a number of electrons per site is 



7TV 
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Appendix D 

Green function for spatially uniform phase 



In this appendix, we explicitly write down the Green function defined by (|3-lf p 
in the spatially uniform case: 



oj (Hkv cos k{£ — y) + A sin k{£ — yj) sin kx 
kv 2 D(u) 

to [Hkv cos k(£ — x) + A sin k(£ — x)) sin ny 
kv 2 D{uj) 

u (Hkv cos kx + A sin kx) sin k{£ — y) 



(x < y) 



(D-l) 



(x > y) 



g—(x,y;u) 



kv 2 D{oj) 

uj(Hkv cos Ky + A sin Ky) sin k{£ — x) 



kv 2 D{oj) 



(x < y) 
(x > y) 



(D-2) 



g+-(x,y,uj) = g- + (y,x;uj) 

hjj 2 s\tlk(£ — y) sin^x 



kv 2 D{uj) 

[Hkv cos Ky + A sin ny) {Hkv cos k(£ — x) + A sin k(£ — x)) 

Hkv 2 D(uj) 



(x < y) 



(D-3) 



(x > y) 



where k = y (Huj) 2 — A / (hv) and D(u) = Hkv cos + A sin 

Appendix E 

Stability of fixed points 



In this appendix, we show that nontrivial solutions of (|4-9p are more stable 
at constant current than those at constant bias voltage. Then, we prove that, at 
constant current, the zero-temperature ordered phase with A given by (I5-I3j) is 
stable. 

Firstly, we note that the stability indices xv a t constant bias voltage and xi a t 
constant current differ by 



-r( 9S 



Xi{A) - xv{A) = -\A — ^= — 



dV 



dJ 



dA 



dJ 



dV 



(E-l) 



As easily seen, we have 
eT 



de 



e sinh ( ^ 



sinh(f)V^y2 J A ^ e2 _ Z 2| cosh 



cosh 



T 



)}' 



> 
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dl\ G /■*"<= Je 2 -A 2 1 + cosh (If) cosh (t. 



cosh (§r) + cos ^ (t) } 



de — , < 



fj ydAjv ja e v^ 2 -^ 2 {^(^) +cosh W} 

Thus, xi {A) > xv (A) which implies that the phase is more stable at constant 
current than at constant bias voltage. 

Now, let us study the stability of the ordered phase given by f)5- 13j) at constant 
current. It is easy to show 

where r = \eV/2A\ > 1. Thus, we obtain the desired result: 

X, . 2A(- - l) + ^-{r - I - r*^HI} 



-Ti 2 J r z -1 J 1 I x 2 x ) 



2A 



frw c \ 2Xr f r , 1 + 3x 2 
2A ; _ / dx—— > 



{hw c y - \a\ 2 

huj c 

(hiu c ) 2 - \A\ 2 * ' r2 ~ x h ^ x 2 (x 2 + l + xVx 2 -1) 

Appendix F 

Ginzburg-Landau expansion coefficients 



In this appendix, we list up the coefficients Ki and K4 introduced in (|5-23[) when 
hw c » T. By letting hio c /T -> 00 in S(A, V, T) - 5(0, V, T) and Taylor-expanding 
the result with respect to A/T, we obtain the desired expansion: 

°-2A + 2 5(0 ' F ' T) + 2 " 2A -^ + ^ ' (F,1) 

where -K2 and K4 are functions of eV/ (2T) defined by 

d ( sinht 

^2 

^4 



o t dt \t(cosh(eV/(2T)) + cosht) 
00 dt d (1 d / sinht 







t cftlt dt\t(cosh(eV/(2T))+cosht)// ' 
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